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ABSTRACT 

Since type la Supernovae (SNe) explode in galaxies, they can, in principle, be used as the same 
tracer of the large-scale structure as their hosts to measure baryon acoustic oscillations (B AOs) . To 
realize this, one must obtain a dense integrated sampling of SNe over a large fraction of the sky, 
which may only be achievable photometrically with future projects such as the Large Synoptic Survey 
Telescope. The advantage of SN BAOs is that SNe have more uniform luminosities and more accurate 
photometric redshifts than galaxies, but the disadvantage is that they are transitory and hard to obtain 
in large number at high redshift. We find that a half-sky photometric SN survey to redshift z — 0.8 
is able to measure the baryon signature in the SN spatial power spectrum. Although dark energy 
constraints from SN BAOs are weak, they can significantly improve the results from SN luminosity 
distances of the same data, and the combination of the two is no longer sensitive to cosmic microwave 
background priors. 

Subject headings: cosmological parameters — distance scale — large-scale structure of universe 



1. INTRODUCTION 

Type la supernovae^ (SNe) have become a ma ture tool 
for studying the cosmic expansion history (e.g ., iPhillipsI 
Il99a iRiess et al.lll998l : iPerlmutter et al.lll999| i. A num- 
ber of SN surveys, such as the S loan Digital Sky Survey 
(SDSS) II (Trieman et al."20078D. the Supernova Legacy 
Survey (Asticr et al. 2006), and t he ESSENCE Super- 
nova Survey (Mik naitis et aL|[2007[ ) are being carried out 
to improve the statistics and our understanding of the 
systematics. Moreover, the SN technique will be an in- 
tegral part of almost every proposed dark energy survey 
inc luding the L arge Synoptic Survey Telescope^ (LSST, 
see lTvson| l2006') and the Joint Dark Energy Mission. 

The conventional SN technique, measuring only the 
relative luminosity distance, is subject to degen- 
eracies between cosmological parameters. For exam- 
ple, the SN constraint on the dark energy equation-of- 
state (EOS, w) parameter Wg., as defined by w{z) — 
wo+Wa,z/{l-\-z), is sen sitive to the prior on the mean cur- 
vature of the universe (|Lindeill2005l : lKnox. Song, fc ZhanI 
[200i ). The reason is that the response of the relative 
distance to a variation in Wa resembles that to a vari- 
ation in the mean curvature and that the S N technique 
lacks the calibration of absolute distances (Zhan fc Knoxl 
I2006D . Even for a flat universe with w{z) = wq, 
the SN constraint on wq can be tightened consider- 
ably if the matter d ensity is known to high precision 
(jFrieman et al.ll2003 ^ Such priors may come from other 
techniques, such as the cosmic microwave background 
(CMB), weak lensing, and baryon acous t ic oscillations 
jBAOs, 'Eisenstein, Hu. & Tegmark' '1998'; 'Coorav et al.! 
2001; Blake & Glazcbrook 2003; Hu & Haiman 2003; 
Lindei1l2003t iSeo fc Eisensteinll2003l) . It has indeed been 
demonstrated that the latter three techni ques are highly 
complementary to the SN Dl technique (jFrieman et al.l 

^ We consider only type la SNe in this Letter. 
^ See ,http : //www . Isst ■ org [ 



[2003tlSeo fc Eisensteinl[20?)3l : iKnox et al.ll200fif ). 

Since SNe explode in galaxies, their distribution bears 
the BAO imprint as well. To measure the SN spatial 
power spectrum, one needs the angular position and red- 
shift of each SN, not its luminosity. Hence, the SN 
BAO technique does not suffer from uncertainties in the 
SN standard candle, which constitute the largest un- 
known in the Z3l measurements. Nevertheless, the nar- 
row range of the SN intrinsic luminosity reduces the effect 
of Malmquist-like biases and luminosity evolution. The 
SN rate traces the mass an d star formation of t he host 
galaxies with a time delay () Sullivan et al .1120061 ). This 
means that SNe have a different clustering bias than 
galaxies that are selected by their luminosity or color. 
Finally, SNe have rich and time- varying spectral features 
for ac curate estimation of photometric r e dshifts ( photo- 
zs) (iPinto. Smith, fc GarnavichI [200l IWa^ 120071 : 
I Wang. Naravan. fc Wood-Vasevll2007| ). which is helpful 
for measuring BAOs from a photometric survey. 

There have been discussion s of using the SN weak 
lensing magnification (jCoorav. Holz. fc HutereJ 
l2006t) and nearby SN pe c uliar velocities 
(jHannestad. Haugboelle. fc Thomsenl l2007| ) to probe 
the large-scale structure. We focus on photometric SN 
surveys for BAOs in this Letter and note in passing 
that the SN weak lensing technique is more limited 
by s h ot noise than the SN BAO technique (| ZhanI 
I2006t IZhan fc Kn"m3 l2006f) and that the SN peculiar 
velocity technique requires precise redshift and distance 
measurements. 

For the BAO technique to be useful, one must survey 
a large volume at a sufficient sampling density as uni- 
formly as possible. Although SN events are rare, the 
spatial density of SNe accumulated over several years 
will be comparable to the densities targeted for future 
spectroscopic galaxy BAO surveys. 

2. PHOTOMETRIC SUPERNOVA SURVEYS 
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We assume two photo-z SN survey models: a shallow 
one (S20k) that covers 20,000 deg^ to z = 0.8 for 10 
years, and a deep one (D2k) that covers 2000 deg^ to 
z = 1.2 for 5 years. The S20k data may be extracted 
as one of the many products from the proposed LSST, 
and a dedicated ground-based SN survey will likely be 
sufficient to produce the D2k data. Since there is no 
precedent for such photo-z SN BAO surveys, we have to 
make crude estimates for the fiducial survey parameters 
and vary them to cover a wider range. 

We calculate the observer-frame S N rate based 
on the rest-frame SN rates (ICappcU aro et al.l Il999l; 
Hardin et all l2000l : 'Pain et al.' 2002; Madgwic k et al. 
2003; Ton rv et al.ir2 003: Blanc et al. 2004; Dahlcn_^LaL 
2004; Nei U et al ll20'0Q) . The resuhing number of SNe per 
steradian per unit redshift per year (observer frame) is 
roughly 



dn 



dQdzdt 



oc e 



3.12z^ 



1 X 



1 z < 0.5 

g-(.-0.5)V2a^ Z>0.5 ' 



where the first term on the r.h.s. fits the observed SN 
rate at z < 0.55, and the additional exponential term 
with a = 0.125 (0.21) cuts off the distribution a,t z ^ 0.8 
(1.2) for S20k (D2k). We take the efficiency with which 
S20k (D2k) will produce well-observed SN light-curves 
suitable for accurate photo-zs from this distribution as 
50% (100%). The accumulated SN surface density is then 
S = 370 and 980 deg-^ for S20k (7.4 miUion SNe in 10 
years) and D2k (2 million SNe in 5 years), respectively. 

We adopt a pedagogical convention and model the 
photo-z error as a Gaussian with rms CTz — crzo(l -I- z) 
and zero bias'^ 6z — 0. We assign for S20k azo = 0.02 
and for D2k azo = .01, which ar e achiev able with simpl e 
photo-2 algorithms (jPinto et al.i r2004: Wa ng et al.ll2007f ). 

Table [T] summarizes the SN surveys. It includes two 
additional quantities: the SN clustering bias b, and the 
cut-off wavenumber /cmax for the BAO analysis, which 
is set to reduce the impact of nonlinear growth (see 
lEisenstein et al]|2007l for recovering high-fc information 
from spectroscopic surveys). 

3. FORECAST METHOD 

We use a modified forecast tool CSWAB'' (|Zhan|[20M) 
to assess the cosmological constraints from the clustering 
of SNe. In summary, the SN power sp ectrum in the ith 
redshift bin Pi{k{) reads ()Seo fc Eisensteinij2003i) 



{l + (3,f,')'b^GMk) (1) 



X exp [-{caz,ik\\/Hif] + Si, 

where the subscript f denotes quantities in a reference 
cosmological model (the same as the fiducial model in the 
forecast), D is the angular diameter distance, H is the 
Hubble parameter, (3 is the redshift distortion parameter, 
H = k\\/k, G is the linear growth factor, V{k) is the 
matter power spectrum at z = 0, and s = is the 
shot noise. The true wavenumbers k±_ and k\\ are related 
to the references by fcj, = kf±Df/D and A;|| = k{\\H/H{. 
Note that for a photo-z survey with cr^ > 0.01(1 + z), the 

^ Any known photo-z bias can be calibrated out in advance, so 
only the uncertainty of the bias matters. 
Available at |http://hzhan. net/soft/. | 



TABLE 1 

Supernova Survey Parameters 





Area 










k 


Survey 


deg2 


z 


h^Mpc"^ 




b 


hMpc-^ 






0.3 


4.2 X 10""' 


0.026 


1.18 


0.17 


S20k 


20000 


0.5 


6.3 X lO-'' 


0.030 


1.30 


0.20 






0.7 


3.4 X 10-"' 


0.034 


1.42 


0.24 






0.3 


4.2 X 10-" 


0.013 


1.18 


0.17 






0.5 


6.6 X 10-"' 


0.015 


1.30 


0.20 


D2k 


2000 


0.7 


7.8 X 10-"' 


0.017 


1.42 


0.24 






0.9 


5.0 X 10-"' 


0.019 


1.54 


0.29 






1.1 


1.8 X 10""' 


0.021 


1.66 


0.34 



Note. — The redshift is central to each bin, and the width 
of each bin is Az = 0.2. 

radial BAO information is essentially lost, e.g., the power 
spectrum (excluding the shot noise) at the fundamental 
mode of BAOs fc|| ~ 27r/150 Mpc^^ is suppressed by a 
factor of 80 or more aX z — 0.6, though a photo-z rms of 
0.003(1 -I- z) would recover that information. 

CSWAB uses the Fisher matrix to estim ate the error 
bounds of the parameters (|Tegmarklll997f l: 



BAO 



a/3 



E 



2 



91nP,(fcf) 91nP,(fcf) dkf 



dqa 



dqp (27r)3 



(2) 



where Vf is the comoving survey volume and {qa} is the 
parameter set, which includes hi, (Jz,i, Szi, and Si of each 
bin and 8 cosmological parameters: wq, w^, (the mat- 
ter density), Wb (the baryon density), 9s (the angular size 
of the sound horizon at the last scattering surface), fik 
(the curvature parameter), Ug (the scalar spectral index), 
and (the normalization of the primordial curvature 
power spectrum). The minimum marginalized error of 

qa is cr{qa) = {F~^)]/a ■ Independent Fisher matrices 
are additive; a prior on qa, ap{qa), can be introduced via 
^aa" = Faa+o'p^{qa)- The fiducial valucs of 5^ , az,i, and 
are listed in Tabled 6zi = 0, and {wq, w^, uJm, Wb, 
Os, l^k, Us, A|) = (-1, 0, 0.127, 0.0223, 0. 596°, 0, 0.951 
2.0 x lO"*^) from the WMAP 3-year results (jSpergel et all 
l2007t l. Unless stated otherwise, we always include fairly 
weak priors ap{lnbi) = 0.3, (Tp(ln A^) = 0.2, (7p(ln0s) = 
o'p(ns) = 0.05, and ap{Szi) = 2~^/'^ap{az,i) — 0.25az,i- 

To show the complementarity between the SN BAO 
and _Dl techniques, we include the Z3l constraints at the 
end of Section [31 The Fisher matrix for SN Z3l is 



1 



Q/3 



n^^z,)^^^^^^^dz. 



Tp^J [51ogZ?L( 



dqa dqp 
z) + M + eiz + e2Z^l p{z\zp)dz, 



(3) 



where tTm = 0.15 is the scatter of the SN apparent mag- 
nitude, the subscript p signifies photo-z space, M is the 
SN absolute magnitude, ei and 62 account for possible 
SN evolution, and p(z|zp) is the probabil ity density of a 
SN at z given its photo-z Zp. Following lAlbrecht et al.l 
(|2006f l. we impose a prior of 0.015 on ei and 62 and let M 
float. We also include z < 0.2 SNe in F^^, because they 
are important for cosmology with Dl (|Lindeill2006( l. 

4. COSMOLOGICAL CONSTRAINTS 

The marginalized 1 a errors of a subset of the cosmolog- 
ical parameters are given in Tabled The two photo-z SN 
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TABLE 2 

Marginalized la Errors on Selected Cosmological 
Parameters from SN BAOs 











lna;m 






Surveys 






Wp 




xlO-2 




S20k 


0.85 


2.9 


0.22 


23 


38 


1.8 


S20k+Planck 


0.66 


1.8 


0.12 


0.83 


0.91 


0.55 


D2k 


0.96 


3.0 


0.19 


25 


42 


2.2 


B2k+Planck 


0.74 


2.0 


0.13 


0.85 


0.91 


0.60 



Note. — The error of the pivot EOS uip equals that of 
a cons tant EOS, and the D ark Energy Task Force figure of 
merit l lAlbrecht et al.ll2006l) equals 0.052[a{w^) Xa{wp)]-'^. 

BAO surveys place rather weak constraints on the dark 
energy EOS even with Planck priors. Nevertheless, they 
provide moderate constraints on the matter density ujm, 
baryon density Wb, and curvature parameter f2k, which 
are helpful for the SN Dj^ technique. The smaller area 
of D2k is compensated by its greater depth and better 
photo-zs, so that D2k performs nearly as well as S20k. 
Spectroscopic BAO surveys of similar characteristics can 
reduce the error on Wm, Wb, and fJk by a factor of 2 and 
significantly more on wq and w^. 

The baryon signature has been detected at ~ 3 cr level 
[and cr(lnti;ni) ~ 10%] from SDSS L uminous Red G alax- 
ies, both spectr oscopically (Eiscnst ein et al.l I2005D and 
photo metrically (jBlake et al., .2007; .Padmanabhan et al.l 
I2007D . These detections assume a flat universe with a 
cosmological constant and a fixed scalar spectral index 
Us- Under the same assumptions, S20k BAO can con- 
strain to 8% and uj^ to 1 5%. If Wb is fixed as well 
(as in lEisenstein et all [20051) . S20k BAO can achieve 

(T(lnWin) = 1-5%. 

Theoretical uncertainties in the redshift distortion, 
nonlinear evolution, and scale-d ependent clustering bias 
can b e imp ortant to BAOs (|Seo fc EisensteinI 120051 : 
I White! 120051 : iGuzik et al.ll2007f) . For a simple test, we 
replace the l inear redshift distortion factor (1-1- 
(|Kaiseij[l987l ) in equation (P) with 1 + 2e/3/x2 + fP'^n'^, 
where e = / = 1 are parameters accountin g for our un- 
certa in knowledge of the redshift distortion (jScoccimarrol 
I2004D . We take priors crp(e) = 0.05 and crp(/) = 0.1 and 
find a less than 1% change to the results in Tabled Spec- 
troscopic (galaxy) surveys with the same redshift distri- 
bution and coverage will see wq and Wg_ errors doubled 
in this test, because they have more information to lose. 

Figure [1] explores the dependence of the SN BAO er- 
ror product <7{wa) X (T(wp) (EP) on the photo-2 rms Cz 
and prior ctp(ctz). At a large CTz, there is little radial 
BAO information, so that the Fisher matrix F^^O scales 
roughly with the fc-space volume, which is proportional 
to HjcUz- This leads to a scaling EP oc when 
and are known acc u rately (das hed li n es), c onsistent 
with iSeo fc EisenstehJ (|2003( ) and IZha^ (|20M) . When 
Cz ^ 0.008(1 -I- z), the radial BAOs become available, 
and, hence, the EP slope steepens at smaller a^- 

Figure[2]shows the EP contours (solid lines) of the shal- 
low survey as the projected SN number density and sur- 
vey area vary. Over-plotted in dotted lines are contours 
of the total number of SNe. If the number of SNe is taken 
as a crude proxy for effort, one can optimize the survey 
(in terms of dark energy constraints) by searching the 
minimum EP on the constant-effort curve. This means 




0.02 0.04 

Fig. 1.— The EP a[wp) X a{w^) of S20k (thick lines) and 
D2k (thin lines) as a function of the rms photo-x error a^. The 
priors on the photo-z biases are taken to be 0.5(T2 (solid lines) and 
0.05(72 (dashed lines), which correspond to calibrations with 4 and 
400 spectra per redshift bin, respectively, in the Gaussian case. To 
reduce the dimensions, we peg the prior on the photo-z rms to that 
on the photo-z bias: a-p(pz) = V2ap(dz). For comparison, LSST 
weak Icnsing, galaxy BAOs, and the two combined will achieve EP 
~ 0.01, 0.02, and 0.002, respectively (iZhanH2009) . 




E (lOOdeg-2) 

Fig. 2. — The EP cr(uip) x (T(uia) (solid lines) as a function of the 
survey area and SN surface density E (accumulated over 10 years). 
The EP contours are spaced logarithmically. The SN distribution 
is scaled from S20k (solid circle, 7.4 million SNe) by S. Photo-z 
parameters are the same as those of S20k, and the Planck priors 
are assumed. The total number of SNe is held constant along each 
dotted line, which is spaced at factors of \/2, 

that for E > 200deg~^ and area less than 30,000deg^, 
one would always choose the maximum survey area pos- 
sible for SN BAO as opposed to accumulating more SNe 
in a smaller area. 

Finally, we illustrate the complementarity between the 
SN BAO and I?l techniques in Figure [3l Although the 
S20k BAO technique (dashed line) does not place use- 
ful constraints on dark energy, its combination (shaded 
area) with the Dl technique (dotted line) reduces the EP 
by a factor of 5.5 over the DL-alone EP. Moreover, the 
combined result (solid line) is not very sensitive to the 
CMB priors, because the BAO technique can provide ad- 
equate constraints on cosmological parameters, such as 
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S20k^ With Planck 

\ " N SN BAO 

SN D, 




■1.4 



1.2 



-0.8 
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Fig. 3. — Marginalized 1 a error contours of the dark energy EOS 
parameters uiq and tUa from S20k SN BAOs (dashed Hne, EP = 0.22 
with Planck), luminosity distances (dotted line, EP = 0.064 with 
Planck), and the two combined with (shaded area, EP = 0.012) 
and without (solid line, EP = 0.018) Planck. 

oJm and r^k, for the technique. 

5. DISCUSSION AND CONCLUSIONS 

We have demonstrated that photo-z SN data can be 
used to measure BAOs and to constrain cosmofogical 
parameters. The BAO constraints on the matter den- 
sity and the baryon density are sensitive to the pri- 
ors on the curvature and the scalar spectral index but 
not to the dark energy parameters. The dark energy 
constraints from the SN BAO technique alone are not 
meaningful. However, a combination of the SN BAO 
and Z?L techniques reduces the EP considerably, and the 



Planck priors are no longer crucial. The SN BAO re- 
sults in Section [3] are also applicable to photo-z galaxy 
BAOs. We note that photo-z errors are a large uncer- 
tainty for photo-z SN cosmology. Although our assump- 
ti on about them is c onservative compared to the results 
in lPinto et al.l ()2004f ). further studies are needed to make 
realistic forecasts. 

Long and non-uniform cadence may result in uneven 
sampling of the SN spatial distribution. Fortunately, 
LSST will be likely to always catch the SNe (especially 
high-z ones) at their maximum owing to its fast sky cov- 
erage and rapid sampling. Furthermore, the effect of the 
cadence on the SN depth can be simulated and deter- 
mined, and methods of correcting for uneven depths in 
galaxy surveys can be applied to the SN data. 

A spectroscopic SN BAO survey will be impractical, 
because one would have to revisit the sky many times 
spectroscopically over thousands of square degrees to 
catch the SNe that occur at different times. However, 
LSST will be able to obtain SNe in the millions over half 
the sky photometrically. This opens a window for ap- 
plying the BAO technique to SNe and achieving more 
robust constraints with photo-z SN data. Since this SN 
BAO analysis requires no additional observations than 
doing the SN Dl analysis alone, it should be a feature 
of all large-area SN cosmology analyses. Moreover, the 
SN data can also help calibrate the host galaxy photo- 
zs and the photo-z; error distribution o f other galaxies 
through the cro s s-correlation m ethod ([Schneider et al.l 
[2001: [ZhanI [20061 : INewmanI [2008') . 

We thank J. Frieman, L. Knox, and M. Wood-Vasey 
for useful conversations. This work was partially sup- 
ported by a UC Davis Academic Federation Innovative 
Developmental Award. 
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